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In modern turbofan engines, variable geometry has been 
incorporated to improve some off-design performance. Most 
control designs ignore this variable geometry and use fuel 
metering as the primary control input. 

This thesis investigates the use of variable geometry to 
control the engine and, thereby, reduce fuel consumption due 
to transients. Additionally, steady-state trim conditions 
are altered to reduce the static fuel consumption. The non- 
linear transient simulation program is used to analyze the 
steady-state operating condition and develop small perturba- 
tion control limitations. Linear models, both large and 
reduced order, are used in analyzing the effect of controllers 
on system response. A computer program was generated to re- 
duce a large order linear model to a usable size for control 
system development. 

This analysis shows the reduced-order regime dependent 
controllers to be viable and to favorably enhance the quest 


for reducing specific fuel consumption in existing engines. 
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I. INTRODUCTION 


High performance, military aircraft are operationally 
tasked to conduct missions throughout their designed per- 
formance envelopes. To accomplish these missions, the air- 
craft must operate efficiently within its flight regime. 

The typical jet engine is designed for operation at a single 
specified altitude and Mach number combination, e.g., 30,000 
feet at M= 0.9. This design point will give inherently 
good cruise performance and is particularly well-suited to 
commercial transport operations. This engine design method 
forces the user into accepting off-design performance in 
Eypeeal Military applications. Engine performance needs to 
be optimized to efficiently conduct all assigned missions, 
@.g., alr combat maneuvering, ground attack, supersonic inter- 
ception and carrier operations. Adaptive digital control 
methods can be incorporated within the engine to accomplish 
this mission dependent optimization. 

Modern control methods allow the designer to develop a 
control system to regulate every parameter of the parent 
system. He has been aided by computer simulations of non- 
linear systems and linear approximations of those systems. 
The desired control system can be developed parallel to or 
ahead of the controlled system. 

The state-of-the-art turbofan engine utilizes variable 


geometry, in such forms as variable incidence inlet guide 


ie 











vanes, adjustable exhaust nozzle area and airflow bleed-off. 
Hydromechanical controls have been used in the past to im- 
prove some elements of engine performance through positioning 
these variable engine components. Electronic controls, such 
as were introduced in the Pratt & Whitney F1l00-PW-100 turbo- 
fan engine, allow improved scheduling of the variable devices, 
again improving the engine performance. Studies have been 
done on the F100 engine to optimize its performance at a 
point other than its design point [1]. 

Originally, the onboard microprocessing capability of an 
aircraft was intended for weapon system improvement and 
integration. Further computing equipment was added as part 
of the engine control to supervise the hydromechanical con- 
trol system and detect faults. The F100 engine introduced an 
electronic control which allowed use of digital control tech- 
niques incorporating the hydromechanical devices as implemen- 
tors rather than controllers. Digital device improvements 
continue to expand the available computing power available, 
making flight control and engine control integration a dis- 
tinct possibility. As more software space becomes available 
to the engine controller, regime adaptive control will be the 
accepted standard in military applications. 

This thesis looks at using modern, multivariate control 
logic in designing flight-regime dependent control for a 
typical military turbofan engine, the Pratt & Whitney F100. 
Specifically, the attempt is made to improve thrust specific 


fuel consumption, i.e., fuel economy, in the high altitude 


2 








cruise environment. Such improvement in fuel economy is 
especially important in the conservation-conscious decades 


ahead. 
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II. TURBOFAN ENGINE PERFORMANCE CONSIDERATIONS 


A. FLIGHT REGIME FACTORS 

The ability of a military aircraft to perform an assigned 
mission is determined by its flight regime. Figure II-1l is 
a typical flight regime for a modern, high-performance 
Misitary aircraft. 

The confines of the flight regime are dictated by numer- 
ous aerodynamic and propulsive factors. A well designed 
aircraft will incorporate an engine having operating boundar- 
les in excess of those of the airframe. This design technique 
allows for a factor of safety in actual operation as well as 
compensate for installation and equipment aging losses. 

At the low speed, low altitude region of the envelope, 
engine performance is judged in the ability to produce 
sufficient thrust for take-off and landing and to provide 
rapid engine acceleration in the missed approach or balked 
landing. The extreme upper region of the flight envelope 
requires the engine to produce maximum thrust while not ex- 
ceeding turbine temperature limitations. 

Figure II-2 is a mission delineated flight regime for the 
Supersonic tactical fighter aircraft. The desired mission 
performance from this chart must be interrelated with the 
aerodynamic and propulsion limitations of Figure II-l. An 
adaptive control system can be designed to schedule optimum 


performance of the airframe and engine throughout the flight 
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regime. Northrup Corporation is involved with a U.S. Air 
Force project to implement an adaptive flight control system 
in an aircraft such as the F-l16. The same type of approach 
can be applied to engine control and integrated into an 


alrcraft control system to optimize the total package. 


B. TURBOFAN ENGINE PERFORMANCE MEASUREMENTS 

The primary factor with which to judge engine performance 
is the engine's ability to convert fuel to thrust efficiently. 
This measurement, thrust specific fuel consumption, SFC, is 
a function of the airflow through the engine, Moe the bypass 
ratio, b, and the fan and core exhaust velocities, Uo! and 


U.': also termed cold and hot stream velocities, 


M 
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ue 


SFC = ((1+f£)U,' + gu," - (1+8)V] (II-1) 


where M 1s the fuel flow, f 1s the fuel-to-air ratio and V 
is the flight velocity. Figure II-3 shows the F100-Pw-100 
engine in cross section with reference labels. 

The engine designer fixes the airflow and bypass ratio of 
the engine. The exhaust velocities are functions of the ex- 


haust temperatures and pressures, which are controllable 


through fuel metering and variable geometry. 
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where y and Y, are the respective ratios of specific heats, 


ac and Tog the respective exhaust total temperatures, n is 


06 
the exhaust nozzle efficiency, and P./Po6 and P./Pog the 
respective ratios of free stream ambient pressure to dis- 
charge total pressure. 

The standard bypass fan will produce a fixed pressure 
ratio and exhaust temperature, producing an exhaust velocity 
in the cold stream in excess of the flight velocity. This 
increases the engine's thrust without altering the hot stream 
through the engine core (ignoring energy lost in driving the 
fan). The advanced, low-bypass fans that reintroduce the 
cold flow into the augmentor, as in the F100, increase the 
augmentor efficiency and during non-augmented operation, 
allow for increasing Boo" through the use of the variable ex- 
haust nozzle. 

Studies show that Thrust Specific Fuel Consumption, SFC, 
can be decreased through increased turbine inlet temperature 
and compressor pressure ratio. Both of these variables are 
designed into the engine at the engine design point. Engine 
control hardware can be designed to maintain the turbine inlet 
temperature, FTIT, and pressure ratio, Rae at the optimum or 


best achievable level during off-design operations. 


C. ENGINE CONTROL DEVICES 
The mechanical devices available for engine control are 
the hydromechanical fuel control, the variable exhaust nozzle 


and variable incidence inlet guide vanes or stators. The 


ro 





first device has the effect of increasing the fuel-to-air 
ratio and temperatures in the turbine and exhaust sections of 
the engine. The combined effect of increasing f and Se 1s 
increased thrust, but does not necessarily improve SFC due to 
the corresponding increase in the fuel flow. The variable 
area exhaust nozzle has been used for the last two decades in 
most high performance military aircraft. This device in- 
creases the hot section exhaust velocity, Cae thus increasing 
thrust and maintaining near ideal pressure in the exhaust 
section of the engine. 

The latest mechanical devices are the variable incidence 
blading concepts. The first use of these devices was to 
achieve near design mass flow rate through the compressor 
during the engine start cycle. These variable stators would 
be rotated to a starting configuration until the engine reached 
a given rotational speed and then be driven to a fixed normal 
operating position. 

The full advantage of variable geometry blading comes in 
the optimizing of off-design performance. The single stage 


pressure ratio 1s found to be 


U z e 


eee P+ nel - “(tan b. + tan a.)] (II-4) 
Pol RT 51 UL 2 Ee 


where n is the stage efficiency, De 1s the rotational velocity 
©: the Dlading, R is the Universal Gas Constant and T51 the 


total temperature at the front face of the stage. At the 
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design condition, the factor (tan b. + tan ay) PoE nimi Zed , 
thus providing the best pressure ratio. The angle b. 1s the 
relative flow direction upon leaving the rotor portion of the 
machine; ay is the angle at which the flow impinges the rotor 
and is a function of the stator angle (see Fig. II-4). Con- 
trolling the stator angle increases the stage pressure ratio 
in the early stages and increases the compressor (or fan) 
performance over a wide range of operations. 

In the FLOO-PW-100 engine, four mechanical devices are 
used, fuel control, variable exhaust nozzle, inlet guide vanes 
with moveable trailing edge and moveable stators in the first 
three compressor stages. The exhaust nozzle maintains nozzle 
area, expansion ratio and boattail drag, simultaneously, near 
optimum. The inlet guide vanes are used ahead of the fan to 
improve inlet distortion tolerance, improve fan efficiency 
and enhance engine acceleration performance. The variable 


Stators improve starting and high Mach number characteristics. 


D. ENGINE MODELING 

It is not economically practical to build an engine and 
then conduct experimentation to determine its performance. 
Modern technology, in both digital and hybrid computers, 
allows the design engineer to build a simulation model of the 
engine and predict the system response well in advance of 
hardware manufacture. Several generic computer simulations 
are in use, these allow non-specific engine performance analy- 
Slis. Additionally, a computer model 1S generated for each 


new development engine by the engine manufacturer. 
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Figure II-4. Compressor Stage Flow 
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Szuch, [2], provides a detailed description of current 
engine modeling techniques. The Pratt & Whitney F100-PW-100 
engine has been modeled with both a real-time hybrid computer 
[3], and a digital transient simulation [4]. Both programs 
have distinct advantages and applications. Although the digi- 
tal program is not real-time, it provides a large volume of 
data and is easily accessed and maintained. A copy of the 
FLOO-PW-100 Transient Engine Simulation Deck (CCD 1103.3.0) 
was obtained from the YF100 Special Project Office of the 
Aeronautical Systems Division of the U.S. Air Force at Wright- 
Patterson Air Force Base, Ohio. This program was modified 
to be compatible with the Naval Postgraduate School IBM 
370/3033 computer system. 

Designing a controller based on the non-linear digital 
Simulation alone is extremely difficult, if not impossible. 
Miller and Hackney [5], developed high order linear approxi- 
mations, at steady state operating points, from the non- 
linear simulation program. The linear system was modeled in 


the classical state variable form 


Pet GU (CieE= 5) 


Ste 
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The 16 system states, listed in Table II-l, were chosen to 
coincide with the states measured in the non-linear simulation. 


These states include 2 rotational speeds, 3 internal total 


ae 





TABEE i r=1 


Engine States 


xl -- Fan speed, (NL), in RPM 
x2 -=- Compressor speed, (N2), in RPM 
x3 -- Compressor discharge pressure, (Pt3), in psia 


x4 -=- Interturbine volume pressure, (Pt4.5), in psia 
x5 -- Augmentor pressure, (Pt/7m), in psia 


x6 -- Fan inside diameter discharge temperature, (Tt2.5h), 
in Rankine 


7 =-—- Duct temperature, (Tt2.5c), in R 

xo = Compressor discharge temperature, (Tt3), in R 

x9 -- Burner exit fast response temperature, (Tt4hi), in R 
x10 -- Burner exit slow response temperature, (Tt4lo), in R 
xll -- Burnerexit total temperature, (Tt4), in R 


x12 -- Fan turbine inlet fast response temperature, (Tt4.5hil), 
in R 


x13 -~ Fan turbine inlet slow response temperature, (Tt4.5lo), 
mm R 


x14 -- Fan turbine exit temperature, (Tt5), in R 
X15 == Duct exit temperature, (Tt6c), in R 


x16 -- Duct exit temperature, (Tt7m), in R 


24 





Pressures and 11 total temperatures. The inputs to the sys- 
tem are the fuel flow, the 3 variable geometry elements pre- 
viously discussed, and the customer bleed air percentage. 
The outputs include the net thrust, the engine total airflow, 
turbine inlet temperature and the stall margins. Addition- 
ally, two fan exit pressure ratios are determined for some 
of the operating points. The inputs and outputs are listed 
in Table II-2. 

The technique used by Miller and Hackney perturbed each 
State, xl, slightly while holding the other states and inputs 
constant. This allows calculation of the deviations caused 
by that state in the other states and outputs. These devia- 
tions then form the F and H matrices of equations (II-5) 
and (II-6). Perturbing each input, ui, in the steady state 
configuration determine the G and D matrices. The sampling 
time was set at 7 milliseconds. 

The linear models obtained are very good approximations 
of the non-linear system but are often too complex and/or do 
not contain the most convenient parameterization to be used 
Piaepmactical @@esign. A reduction of model order and, per- 
haps, augmentation of the linear model is generally required 


prior to using the model for control design. 


B, GOMPROL DEVICE EFFECTIVENESS 
Before designing an engine control, one must determine 
the effect each input has on the system and the responsive- 


ness of each input. The five control inputs listed in Table 
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u3 


u4 
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y2 


v3 


y4 


v5 


y6 


Tae Li=2 


Engine Inputs 


Main burner fuel flow, lb/hr 


Nozzle jet area, f£t**2 


Inlet guide vane position, deg 


High compressor variable vane position, deg 


Customer bleed flow, % 


Engine Outputs 


Net thrust, lb 


Total engine airflow, lb/sec 


Turbine inlet temperature, R 


Fan stall margin 


Compressor stall margin 


Fan exit delta p ratio (test data) 


Fan emit delta p ratio (theory) 
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II-2 are used in controlling the F100 engine. The customer 
bleed air percentage, u5, is essentially constant in steady 
State operations and can therefore be ignored in the control 
development. A study could be made to determine the merit 
in extracting additional bleed air from the engine as a pri- 
Mary control, and dumping that portion not required by the 
customer systems. 

The intent of this control development is to minimize the 
specific fuel consumption of the engine. At the same time, 
only small variations in turbine inlet temperature, y3, and 
total engine airflow, y2, are desirable and no change in 
thrust is allowable. One method might employ using exhaust 
nozzle area, u2, and fan inlet guide vanes, u3, inputs to 
minimize SFC at the trim condition and then compensate for 
thrust loss by using fuel to balance the equation. 

Figure II-5 plots the change in specific fuel consumption 
due to incremental changes shown in Table II-3 in the inputs 
ul, u2, and u3; input u4, the compressor variable guide vanes, 
exhibit no influence on SFC. Proportional derivatives, l.e., 
slopes were determined by making a linear approximation to 
the curves. It is clear in the figure that nozzle area and 
inlet guide vane inputs (u2 and u3) have the most substantial 
effect on SFC. 

Figure II-6 plots the net engine thrust, yl, versus the 
control input changes. Again u2 has the greatest effect and 
ul and u3 contribute about equally to thrust variations. It 


1s also noted that a linear combination of ul and u3 can be 
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PASE II-3 


Input Perturbations 


The following perturbations are used throughout this 
investigation. They represent the largest expected value, 


Of Physical limitation placed on each control input. 
Fuel Flow, ul + 200heL pyaar 


Nozzle Area, u2 H0..2 se. ft. 


I+ 


Inlet Guide Vanes, u3 4 degrees 


1+ 


Compressor Vanes, u4 4 degrees 


Bleed Air, u5 0.0 percent 


used to offset the thrust change caused by a change in u2. 
This is important in meeting the zero net thrust change de- 
sign requirement. Figure II-7 shows the effect the inputs 
have on turbine inlet temperature, y3. Again, u2 is the 
dominant control input. These figures clearly demonstrate 
the need for minimizing the exhaust pressure mismatch with 
ambient conditions. 

Table II-4 is a summary of the proportional derivatives 
extracted from the preceding figures. Using these derivatives 


the following relationships are found. 


d (SFC) = 2.27E-05 d(ul) + .0655d(u2) 
- .00092d(u3) (II-7) 
d(yl) = woo omental) = 91351. /dataz) 
+ 18.6275d(u3) (II-8) 


ou 





TABLE II-4 


Summary of Derivatives 


Input SFC iinwst (yi) Bey.) 
Fuel Flow, ul 2.2725 mo eeO ze 6 
Nozzle Area, u2 £0655 =155 1... / 235.085 
Inlet Guide Vanes, u3 =O Wa ono 2 75 -2.1641 


U4 was found to have no effect on the parameters of 
interest. 


Mees) = .1l2iz6di(ul) + 283.085d(u2) - 2.1641d(u3) 29) 


If a decrease in u2 of -.2 square feet and an increase of +4 
degrees in u3 are applied, a change of -.017 occurs in SFC. 

The fuel flow, ul, 1s used to balance the thrust equation. 

To accomplish this the fuel is reduced by 370 pounds per hour, 
further decreasing specific fuel consumption. The total effect 
of these three inputs at the trim condition is a 3.5% reduc- 
tion in SFC, no thrust change and a 110 degree decrease in 

the turbine inlet temperature. 

By requiring that thrust level not change, the SFC per- 
formance 1S compared under identical operating conditions. 
Each input can be studied separately to determine its effect 
on the outputs and on SFC. To remain in the small perturba- 
tion regime inputs should not be varied more than 10% of their 
steady-state values. Combinations of inputs are then chosen 
to effect the most desirable improvement in SFC and minimize 


the effect on engine performance. 
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This technique has shown improvement in the specific fuel 
consumption at the operating point chosen in the studies of 
Reference 1. The use of the existing variable geometry to 
improve the fuel efficiency at the high altitude cruise 
operating point chosen for this investigation will be dis- 


cussed ina later section. 
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Pei CON PROu DEVELOPMENT 


A. LINEAR OPTIMAL CONTROL THEORY 


Given the linear system, 
X = FX + GU (III-1) 


consider a linear control law defined by 


U = ~-KX + Up Citi 2 ) 
The original system can be modified to 
X = (F - GK)X + GU, (Titi =—3) 


The response of the modified system, as well as its eigen- 
values, is determined by the effect of the gain matrix, kK. 
For a low order system, the control gains could be found by 
trial-and-error methods and applied to the system to determine 
if the response was acceptable. This method is impractical 
for large order systems. 

Optimal control is achieved by minimizing a cost function, 
J, which is defined as 

i iD st 


mes 2X) (Ue ede (III-4) 


where R) and Ro are symmetric weighing matrices; Ry 1s non- 
negative definite and Ro LS) posttElve @eeinite. “A third portion 


can be added to (III-4) to provide terminal state weighing, 
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ah 
Xx (t,)P,xX(t,) (ii) 


where Py 1S a non-negative definite symmetric matrix. 
It is shown by Kwakernaak and Sivan [8], that the optimal 


maiDut, U,(t), based on the cost function is given by 


Un(t) = (R,)7°G™P(t) (ch) 
thus the knowledge of p(t) solves the regulator problem. A 
2n linear system is formed with the optimal system behavior, 


X(t), and the adjoint variable, p(t). 


E il 
-l? 
XQ (t) F -GR, G XQ (t) 
= (III-7) 
p(t) = Te p(t) 
| | i t 
where 
p(t) = P(t) X) (t) (III-8) 
with 
ae R57 (t)G" (t) P(t) (III-9) 


P(t) is a symmetric non-negative definite matrix that 


satisfies the Riccatil equation 


35 








P(t) = R(t) - P(t)G(t)R>~(t)G" (t) P(t) 


eee (Ee) P(t) coin) 


with the terminal condition of 


P(t,) = Po (III-L1l) 
Under these conditions, as the terminal time, tl, approaches 
infinity, the control law asymptotically approaches a unique, 
stable steady state condition. Thus, the matrix Riccatil 
equation, (III-10), can be rewritten 

ter Z 


0 = R, ~ PGR, Gor + Pre foe (TTT=12) 


RK = @8R>“GP (Ti T= 13} 


Many numerical methods are available to solve the matrix 
Riccatli equation for the steady-state control gains. Some 

of these are covered by Kwakernaak and Sivan [6]. This 

thesis employed a computer solution, OPTSYS4, developed by 
Hall and Bryson [7], at Stanford University. A representative 


output is included in Appendix B. 


ae) CONTROL SYSTEM MODELS 
At this point, an analysis of control models is made to 


establish the form of the controller to be applied. A 
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comparison of small and large perturbation controls must be 


made. Take the linear system 


Sd 
i 


a tO (III-14) 


K 
It 


ent a) Crit — 15) 


This system can be made closed-loop by incorporation of a 


Gentrol law 
U = =- KX . Gi e— 16) 


K can be a time-varying matrix or constant as required by 


the system. A more general control law would be 


els (LEir=l7) 


eos 2 CL (Xx) + J CACY =, 


where A is the diagonal matrix whose diagonal terms are l 
or 0 determining which outputs are used in formulating the 
TaD wNe . Cy. and Se are again matrices associated with the 


states and the outputs, respectively. Additionally define 
u = U, + Rit Cnr L= hs) 
Xx Ne Rt (TIT =59)) 


the 0 subscript indicates trim position/value; RR and R, are 
Functions of u and x associated with time. For small pertur- 


bations the time varying terms, RX and Ri can be ignored and 


ey 






a 


equation (III-17) can be rewritten 


6U = C. 6X + {C_AS6Y dt (III-20) 
x a 


Dropping the delta terminology and taking the derivative of 


this and using equation (III-16), yields 


=e Cex tC AY Giir=21) 
* Sf 
combining terms 


(-K-C) X = C AY (ia 2 2 ) 


Substituting, using the system description of equations 


(III-14) and (III-15) one obtains 
(-K-C_) (FX+GU) = a (TPp—2 3} 


and solving for Cyr 


GC = -K-C,A(H-DK) (F-GK) ~* 


a Ci 24)) 


and using all outputs (A = I) and equally weighing them 


(Cy = I) one has 


es -~K-(H-DK) (F-GK) ~> (ITI-25) 
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However, if the outputs are ignored (A = 0), 1.@e., consider- 
ing the case where small changes are made in the output 


vector, then 
c = - K (III-26) 


aS in the simplified control law of equation (III-16). 

This approach does not apply to the large perturbation 
case. In large perturbation analysis, the time varying 
terms cannot be ignored. The mean values of u and x are as 
in equations (III-18) and (III-19). RY and Ry are the time 
functions that describe the path that the mean values follow 
over the period of the large perturbation. The control law 


becomes 
Sees eee GXi= Rot) J CASY dt (III-27) 


Again, if all outputs are used and equally weighed, then 


taking the derivative 
ce = oo 5 — CoRR. —-C.R t +.C AY (III-28) 
u u x Scat > y 
and combining terms 


~C (F-GK)X +R C, +R Clt = K(F-GK)X + Ry 


+ Rt + (H-DK) X fie laeo) 


So 











If we set 


RUC antl CIEE S18) 


and 


Re = R (IEI-31) 


in effect, assuming a linear relationship between the states 
and inputs. Then 


-C. = -K - (H-DK) (F-GK) 7) 


(III-32) 
aS in equation (III-25) for the small perturbation case. 

An adaptive control system removes the need to model the 
large perturbation response. Linear models can be developed, 
as described in Section II.E, to define the large perturba- 
tion case as a Series of linear relations. The onboard digi- 
tal computer then schedules the control gains based on the 
linear models as the flight conditions pass through the model 
transition points. The rapidity of response inherent to the 
digital computer removes the requirement for modeling the 
total dynamic system, rather, only the dominant response model 
for the selected operating point. The reduced-order linear 
model is far easier to handle in formulating the control sys- 


tem and is derived to match the dominant system characteristics. 
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C. REDUCTION OF LINEAR MODELS 

DeHoff and Hall [8], present a method of reducing the 
order of the linear system based on the dominance of the 
states and the crosscoupling of the states as observed through 
eigenvector analysis. The required order of reduction is 
determined by the number of dominant states appearing in the 
desired control bandwidth. For the F100 study, DeHoff and 
Hall chose a control bandwidth of 0 to 10 rad/sec, corres- 
ponding to the primary control device (fuel flow actuator). 
Their operating point (static, sea-level, intermediate power) 
dictated that the control be designed to modulate thrust, 
providing maximum thrust for take-off performance. At another 
operating point a different parameter may dictate the control 
design and, thus, the primary control actuator and control 
bandwidth. 

The method of decomposition involves finding the eigen- 


vVeeter transform matrix, T, such that 


XT = TA Ciigie— 3 3) 


where A is the diagonal eigenvalue matrix. Then, by defining 


an alternate state vector, Z, where 
x = TZ (ITI-34) 
the original system may be written in modal coordinates as 


2 = iz + tley (III-35) 
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If a control bandwidth is determined, 


the eigenvalues 


falling in that range are the only ones that need to be in- 


cluded in the initial reduced order model. 


The matrix may 


be reordered by constructing a matrix Q, which has al in 


the position corresponding to the eigenvalue to be included 


(column) and the reordered position (row) 


of that eigenvalue. 


Thus, if eigenvalue 5 of the original A matrix was to be re- 


ordered to the first position, a 1 would be placed in the 


(1,5) position of the Q matrix. Now, the reordered matrix 


is incorporated in the linear model 


oz = A'z + gt lau 
letting 
7 eZ 
Zz = j'q tz: + or tpu 
and 
oS ol 


xy metruxe R can be found such that 
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(Eit= 36) 


(IITI-37) 


(121-38) 


(ilps) 


(hii 20) 


(III-41) 





Mhis X' vector “is the reordered State vector, 


elements contained in the subvector X 


i: 


where the 


those states 


associated with the eigenvalues in the control bandwidth 


and X. 


contains all other states. 


The linear system may be rewritten 


Now define 


xt 


3 


RFR “x' + RGU 
HR tx! + DU 
RFR + 
RG 
HR 
t 
aah 1 Pi2 
ae. ---=- 
I 
Por t Foo 
Gy 
G5 
ee 
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(i= 4 2 ) 


(LIIT-43) 


(III-44) 


(ITI-45) 


(III-46) 


(ILTI-47) 


(III-48) 


Cin 2) 





If the states in X) 


then Xo will be essentially zero and 


Xx = Pii*l at Pixs ate GiU 


0 = Poi) ate Po2%5 te G4U 


Solving for X, in equation III-5l, one has 


2 


_ et =i 
Xo = ‘cpm 


This is incorporated into the first equation 


1 


X, = vl, 


l ee oe PP NE SEF 5) Y 
and the output equation becomes 
ee ie Sh ae er.) eee (DEH. Gu 
noe 22 6 Zee 2 


The reduced order model can be written as 


x = Frx 


1 l + Gru 
1 a= Hox, + Dru 
where 
Fr = F ee a 
nea leg 22 “Wi 
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truly model the full order systen, 


(III=-50) 


(III-51) 


(TIiT=-32) 


Ciro) 


Cie ) 


Clr) 


(ee 6) 


(III-57) 








Cmemmcre- £5255 G5 Chie 5 6) 
eee: & mer (III-59) 
qm DD): 
Dr = ov-H.F.21c (III-60) 
pe? 2 


This technique was used to develop the program REDUCR2, 
documented in Appendix A. REDUCR2 was used to develop the 
reduced order models used in the control development portion 
of this thesis. Sample input and output of REDUCR2 are 


listed in Appendix A. 


D. REDUCED ORDER MODELS 

Three reduced order models were generated with the REDUCR2 
computer program (Appendix A). The criteria for selecting 
the desired states was that used by DeHoff and Hall, a con- 
trol bandwidth of 0-10 Hertz. The eigenvalues of the full 
loth order system were found using the OPTSYS4 computer pro- 
gram and then associated to the system states using modal 
analysis. Table III-l gives the results of this modal analy- 
Sls. The five states associated with the control bandwidth 
are fan speed (xl), compressor speed (x2), augmentor pressure 
(x5), burner exit temperature (x10), and fan turbine inlet 
temperature (x13). 

A fifth order model was generated using these states. It 
was noted that only three of these states (the two rotational 


speeds and the pressure) were dynamically significant, i.e., 
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TABEEB IITI-1 


Operating Point Eigenvalues 


epate Bovico cal Property Bag Chwawe 

x4 PO04.5 ~540.23 

X3 P03 =156.432 

ee T04 -49.69 + j13.99 

X16 EOF 

X9 TO4Hi - 54.12 

al 2 TO04.5Hi - 42.61 

x7 moe oC - 22.59 

X14 TOS alee) sf 2c Ok 

X15 TO6C 

X6 Oe, oH Scenic 

X8 03 - 14.49 

X5 P07 Sian) 45) 

XL Nl es Od Comero lL 
band- 

x3 T04.5L0 - 1.98 width 

X2 N2 Sod 

X10 NO4LO - ~647 
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XL 


X2 


X5 


X10 


os] 


X13 


X16 


vi 


a 


Model #1 


Sao |. 


Zoe O 


oO 


94.48 


)1..034 


Tapes LLL =2 


Model #3 


er 9 


PAO OO 


eo 


boa oe 


23.04 


174.8 


Die oe 


Reduced Model Comparisons 


Model #2 


Oo 


26.0 6 


soe) 


Sie0 


mall 


Se 


Pas 


33% 


Zee 


156 


Sse 


Order 


41 


OZ 


“o26 


.86 


48 


oS 


03 


299 
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represent energy storage in the form of torque and pressure. 
A third order model was generated using these three states 
and a comparison made. In this initial analysis these two 
models produced virtually identical response with respect 
to the key outputs thrust (yl) and turbine inlet temperature 
(y3), see Table III-2. The OPTSYS4 program was used to deter- 
Mine the transfer functions associated with these two outputs 
for each model. It was found that the fifth order model con- 
tained zeros near the poles at -1.960 and -.65l, the eigen- 
values associated with the two temperatures in the fifth 
order model. This caused these two poles to have no effect 
on the response of the model. Thus the third order system 
sufficiently models the control bandwidth. 

In both models, the turbine inlet temperature (y3) response 
did not closely model that of the full order system. There- 
fore, a third model was formed using the three states of 
the third order model plus the conjugate pair of poles repre- 
senting the turbine inlet temperature (xll) and exit tempera- 
ture (x16). This model gave excellent thrust and temperature 
response (yl and y3, Table III-2) due mainly to the inclusion 
of the turbine inlet temperature as a system state. Table 
III-3 lists the three models, their states and associated 
eigenvalues. 

Although the conjugate pair of eigenvalues lies outside 
the control bandwidth, their inclusion markedly improves the 
output correspondence with respect to the full order system, 


see Table III-~-2. Further analysis will use the third order 
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Model #1 


XL 


States 


X5 


X10 


X13 


ele) 2 


Eigen- 

-3.392 
Values 

-1.418 


ete 


= 969 


TaBes Lit—3 


Reduced Model Eigenvalues 


Model #3 Model #2 

X1 X1 

X2 X2 

eS X5 

el 

X16 

a5. 65. = -6.947 
17 a4 -3.282 

- 1.388 -1.389 

- 6.909 

Sa). cae 
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model (model #2) and the fifth order augmented model (model 


#3). These models are listed in Tables III-4 and III-5. 


Pee VoRLELICATION OF REDUCED ORDER MODELS 

The reduced order model must be sufficiently representa- 
tive of the original system for the designed control to pro- 
vide the desired full system response. Weinberg and Adams 
[9], used the F100 non-linear simulation to develop a linear 
17th order system at the static, sea-level point at partial 
power settings. Their reductions to 5th and 3rd order and 
Similar reduction to 5th order by DeHoff [10], from Miller 
and Hackney's loth order model show excellent correlation to 
the non-linear simulation. 

The operating point chosen for consideration in this 
thesis is 30,000 feet, Mach = 0.9 and power level angle of 
67 degrees. This corresponds to a near optimum cruise point 
and is a condition frequently encountered in routine opera- 
tions. Section III.D lists the states and associated eligen- 
vaiues at this operating point and discussed several reduced 
order models. The dominant elgenvalues were considered to 
be those which affected the models fast response and had 
dynamic significance to the engine's physical characteristics. 
Again, this implies a control bandwidth of 0-10 rad/sec, how- 
ever, other models, using eigenvalues outside the control 
bandwidth, were considered in determining the best reduced 
model for deriving the control laws. 

One method to analyze model correspondence is to use the 


Beeemploreacross the control bandwidth. Figure IITI-1l is a 
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Bode comparison of the DeHoff reduced model and the Miller 
and Hackney loth order model at static, sea-level, inter- 
mediate power for the 0 to 10 rad/sec control bandwidth. 

It is clear that through the 0-10 Hertz bandwidth the two 
models show very close correspondence; this similarity breaks 
down in the higher frequency region, but the response time of 
the eigenvalues outside the control bandwidth is considered 
to be too short to affect transient behavior. The phase 


response, not shown, follows the same trends. 
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IV. ANALYSIS OF CONTROL SYSTEMS 
A. COST FUNCTION WEIGHING FACTORS 
The control law for the engine regulating system, 


tee UL + C(X-x) + f CJA(Y-¥q) dt (IV-1) 


0) 


can be reduced as shown in Section III.B to 
U = Uj - KX (IV-2) 


when considering small perturbations and ignoring the out- 

puts from the system. In most circumstances, the Ua! the 

trim condition, would be considered as zero for steady state 
analysis. The linear models developed by Miller and Hackney 
[5], set all inputs to zero at the steady-state operating point 
for which the model is defined. 

The K matrix 1s the solution to the optimum regulator 
problem discussed in Section III.A. The OPTSYS4 computer 
solution used in this thesis uses the cost function 

J = 5 f (vTayy - UR,Uldt (IV-3) 
where oy is the weighting matrix for the output vector and 
replaces the need for the Ry MaGCrix in@éqW@ation (III-7). 

Kirk [11], discusses the effect that weighting matrices 


have on system performance. In general, the larger the 
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magnitude of the matrix element the faster its associated 
parameter (system state, input or output) stabilizes. An 
additional factor in determining the magnitude of the diagonal 
matrix elements is the magnitude of their associated param- 


eters. Given the cost function 


ee AG = 
De ae cle (IV-4) 
where 
ale 
my 0 sel) 
a 
Ry = : x = (IV-5) 
0 = ev 
b 


sos & (HS)? + (S47) at (IV-6) 


The desire is to minimize the cost function, J, with 
respect to time. Thus, taking the time derivative and 
setting it equal to zero 


+ (—=) = Q (iV= J) 


This is the solution which produces a minimum cost. The 


values of a and b are the largest expected values of xl and 
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x2, respectively, and represent an equal effect on the cost 
function by each state. Should either a or b be decreased, 
the optimal control developed would have to increase the 
control input to the other state to minimize the cost func- 
tion. A nominal fluctuation limit of +5% will be considered 
as a small perturbation in the initial control analysis. 
Designing a control to minimize the thrust specific fuel 
consumption (SFC) of the engine requires limiting the thrust 
and the fuel flow fluctuations. The weighting matrix ele- 
ments corresponding to these two parameters are, therefore, 
to be increased with respect to the other parameters. Looking 
at only the small perturbation response, it 1s noted that 
very little change in the total engine airflow (y3) 1s ex- 
pected and, as before, there 1s no change in customer bleed 
alr requirements (u5). These two parameters can be neglected, 
as their response characteristics have no effect on the desired 
SFC response. To demonstrate the effect of the different ele- 
ment magnitudes on the control gain matrix, the OPTSYS4 pro- 
gram was run for the 3 cases delineated in Table IV-l. Case 
I was developed to allow each output and input, with the excep- 
tion of the two pressure ratios (y6 and y7) and the bleed air 
percentage (u5), to effect equal contribution to the cost 
function. The diagonal elements of both AY, and Ra normalize 
the squared value to the order of one. Case II uses all of 
the same elements, except the thrust (yl) and the fuel flow 


(ul) are weighted more heavily, to the magnitude of ten. This 


a7 





ee 


tape. Ly =] 


Case Definition 


A and R. are diagonal matrices, only the significant 


elements are listed. 


Case I 
AY, = CeO One eel: -0001L1 100 100 le-8 1le-8) 
Ro = (20004 100 Onli -OLl le-8) 

Case II 
AY, = (.OLl 1 .0001L1 100 100 le-8 le-8) 
Ro = (.OL 100 On: -Ol1 le-8) 

Case IIl 
AL, = (1 1 .0001 100 100 le-8 Ile-8) 
R, = (1. 010 SO. -O1 le-8) 


heavier weighting should reduce the fuel input to the system 
and moderate the thrust fluctuations as the control returns 
the system to the steady state condition. Case III increases 
the weighting of yl and ul to achieve a normalized magnitude 
of 100. This extremely heavy weighting should eliminate fuel 
from the control and reduce thrust variations to absolute 
minimum. The whole intention of the controller design is to 
damp out state perturbations without using fuel or altering 
the thrust level, i.e., making no change in specific fuel 


SeMmctmetLton to GEEect control. 
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These three cases were implemented on the two reduced 
order models and simulations run for a perturbation of -50 
RPM in fan speed (xl). The fifth order model returned the 
fan speed to 95% of steady state in .5 seconds using Case I, 
but increased fuel flow by 10 pounds per hour initially and 
required an increase in fuel consumption of .2% over the 
half second period. For Case II, the fifth order model re- 
turned to 95% of the steady state fan speed in approximately 
.3 seconds, while using almost 80% less fuel than the Case I 
controller. This 80% reduction in fuel consumption indicates 
that the increased weighting of Case II has accomplished 
exactly what was desired. Case III was run for the fifth 
order model and, although virtually no fuel was used in the 
control, the system experienced oscillations in fan speed of 
+8 RPM through the 3 second simulation run. Additionally, 
the Case III controller caused rapid changes of +4 degrees in 
the compressor guide vanes (u4), the regulatory limit imposed 
for the simulation, and also large fluctuations in the inlet 
guide vanes (u3). One would achieve results of similar magni- 
tude for each of the cases assuming a +1% perturbation of all 
states in the reduced-order model. Case I is too fuel de- 
pendent and causes large changes in u3, Case III also requires 
large changes in U3 that could cause a slew rate request larger 
than that available from the control actuator. 

In modeling the control, the actuator and physical limi- 


tations of each control must be considered and appropriate 


Do 





limitations placed on the control law to avoid overshoots. 
Table IV-2 gives the actuator limits and maximum rates 
applicable to the F100 engine [4]. A stepping limit of .0Ol 
seconds is assumed in determining the rate limits and fuel 
flow fluctuations are limited to +10% of the steady-state 
value. No compensation is allowed for hysteresis or non- 
linear dynamics of the actuators; these are considered 
Minimal in the small perturbation case, but would need to 


be incorporated for large perturbation control modeling. 


TABLE IV-2 


ARetuator Limitatzvons 


iapuat Maximum Minimum Rae lane 

Fuel flow VEsoo- lb nr 450 lb/hr Mae 00 Lb7hr/sec 
Nozzle area G.4 SQ ft. Zeomog. ft. S20 SQ-ft./sec 
Inlet guide vanes 0 deg. -40 deg. 48 deg/sec 
Compressor vanes 4 deg. -40 deg. 40 deg/sec 


It is evident that the Case III weighting is not accepta- 
ble because of the oscillations and the rapid, large amplitude 
Fluctuations in the guide vane positions (u3 and u4). Case I 
is also inappropriate for the specific fuel consumption mini- 
mMLZing control since it uses fuel as the primary control in- 
put. For these reasons, Case II is chosen as the weighting 
for all control modeling to be implemented in later sections. 

Once the integral portion of equation (IV-1) and the feed- 


back gains are determined, the constant portion, on must be 


6 0 





analyzed. In most situations this term would be neglected. 
The steady state analysis is now performed to determine if 
the specific fuel consumption can be reduced by making small 
changes to the trim position of the variable geometry. An 
analysis of control changes versus SFC, as was done at the 
static sea-level condition in Section II.D, must again be 
conducted. A minimum SFC obtainable relates to certain input 
values in U., which is incorporated into the control law, im- 
proving the steady-state performance and controlling the 
fluctuations of the system which will have the greatest effect 
Saacne orc. 

As was done at the static sea-level operating point, a 
small perturbation in each input was made and its effect on 
the SFC, thrust and other engine conditions tabulated. The 
goal of these analyses is to develop a trim input, Us that 
improves the fuel economy at steady state without causing 
Significant change in the operating conditions (both states 
and outputs). 

The effects of the control inputs on several parameters 
are plotted on Figures IV-l through IV-4. Figures IV-l 
shows, again, the large effect nozzle area (u2) and fuel flow 
(ul) have on specific fuel consumption, the effect of inlet 
guide vanes (u3) is drastically reduced from that observed 
at static, sea level operation. Figure IV-2 indicates the 
offsetting effect of ul and u2 and again indicates a zero 


net thrust control can be found. Fan speed (xl) and compressor 
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speed (x2) effects are shown in Figures IV-3 and IV-4, respec- 
tively. At this operating point the inlet guide vane posi- 
tion (u3) may be most useful in controlling the fan speed 
fluctuation and thereby stabilizing the fan stability margin 
(v4). A sample control would close u2 by .2 Square feet and 
reduce fuel flow by 57 pounds per hour. This combination 
produces a 1.4% decrease in specific fuel consumption and 
results in a 4% reduction in fan speed. A trim condition 


Gentrol input, Ua! would be 


B. REDUCED ORDER SIMULATION ANALYSIS 

The model developed in Section III.C and Section III.D for 
the high altitude flight operating point was used to initially 
analyze the control law developed in Section IV.A. A reduced 
order controller is simpler and cheaper to implement through 
a digital computer. It requires less computer space and far 
less computing time than the full-order linear model developed 
control ler . 

A comparison of the three cases for which gain matrices 
were found indicated that Case 2 gave the most desirable re- 
sults. The Case 2 weighting matrix was applied to the third 


order reduced system and perturbations applied to various 
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states. The same weighting matrix was applied to the fifth 
order model (#3) and similar perturbations applied. Figure 
IV-5 shows the system response to a 50 RPM perturbation of 
the fan speed. The fifth order system has a settling time 
of approximately .45 seconds whereas the third order system 
has one of 1.5 seconds. This variation is accounted for by 
the two fast response poles in the fifth order model, 
(5 ,37.14). 

These results must be compared to the l6th order system 
to determine which model most closely approximates the full 


order system response. 


C. FULL ORDER SIMULATION ANALYSIS 

The proof of the viability of a reduced order controller 
1s its ability to control the full order system in the same 
manner as predicted in the reduced order model. The two 
reduced models analyzed in Section IV.B were applied to the 
sixteenth order linear model. Additionally, the Case 2 weight- 
ing matrix was applied and a full order controller was found 
and implemented. All simulation was done using the IBM Con- 
tinuous System Modelling Program, Model III. 

All three models were compared with similar perturbations 
applied. Figure IV-6 compares the full order and fifth order 
response to the 50 RPM speed perturbation. Figure IV-7 com- 
pares the full order and third order response to the same 
perturbation. Figure IV-8 compares the two reduced order 
controliers and the full order controller for the same fan 


speed perturbation. 
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The results show satisfactory agreement between the re- 
duced order and full order implementation of the reduced 
order controller. This proves the viability of the con- 
troller and implies that comparable results will be found 
if the control is implemented in the non-linear dynamic 
system. The fifth order controller, with two poles outside 
the control bandwidth, most closely models the response of 
the full order controller. As was seen previously, the 
inclusion of the conjugate pair of poles dramatically im- 
proves the response time of the model and provides a much 
closer approximation to the full order system response. 

This indicates that the fifth order model is sufficient 
for the control development process at this operating point. 
The two fast response poles outside the control bandwidth 


must be included to achieve this close simulation. 


V2 





V. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

The study of implementing a reduced order controller 
that is flight regime dependent shows that at a specific 
operating point an improvement in specific fuel consumption 
can be made. The conclusions reached at the end of this 
research are: 

1. The variable geometry existing in the state-of-the- 
art turbofan engines can be used to further improve the 
specific fuel consumption; 

2. A reduced order controller can be found and imple- 
mented that concentrates control authority on the variable 
geometry rather than the fuel control, thus furthering the 
fuel savings; and 

3. A similar controller can be found and implemented at 
any given operating point throughout the flight regime, these 
controllers could then be organized into a single digital 
machine to provide a regime-dependent engine control to im- 


prove the airframe/engine interaction in all flight conditions. 


B. RECOMMENDATIONS FOR FUTURE STUDY 

One obvious extension of this work is to find further 
reduced order controllers. 

Additionally, work can be done to fully implement these 
reduced order controllers in the non-linear model using 
either the Pratt & Whitney transient simulation digital model 


or the hybrid computer model. 


me 








The effect on fatigue life and life cycle cost due to 
the "fuel saving" control schemes presented was not deter- 
mined. It is therefore recommended that these factors be 
explored in order to fully analyze the usefulness of these 


control schemes. 
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